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I. INTRODUCTION 

A clean extraction of excited hadron masses from a 
lattice QCD simulation is a serious challenge. However, 
an ab-initio determination of properties for excited light 
hadrons would provide highly interesting information on 
the chiral dynamics of QCD. 

Excited hadrons are rather non-trivial objects to study 
on the lattice. One of the main reasons for difficul- 
ties is the fact that excited states appear only as sub- 
leading contributions in Euclidean two-point functions. 
Although a variety of other approaches has been tried, 
including Bayesian methods [3, [3, Q, an NMR-inspired 
blackbox method 0, Q and evolutionary fitting tech- 
niques the most powerful method is probably the 
variational approach 0, HJ • The reason for its power is 
the fact that in the variational method not only a single 
correlator is studied, but a whole matrix of correlation 
functions. Consequently more information is extracted 
from the system. 

The successful implementation of the variational 
method hinges crucially on the set of basis interpolators 
that are used in the correlation matrix. A particularly 
important criterion is that the interpolators have a large 
overlap with the physical states in a given channel, i.e., 
with both, ground- and excited states. In this article 
we build on earlier work 0, [l^, HH, [Hj]> where Jacobi- 
smeared quark sources with different width were used to 
construct hadron interpolators that allow for nodes in 
their radial wave function. To construct a richer set of 
interpolators, we now also include derivative sources for 
light-quark spectroscopy. 

Interpolators with derivatives have been wi dely used 
for heavy quark systems (see for example (15 , [lj, [TEl]). 



They have also been applied by Burch et al. [H, [13] for 
light mesons and by Lacock et al. [H[ to the study of both 
orbital excitations of mesons and hybrids. Note that the 
approach in [l8| differs from ours. There interpolators 
are built with quarks displaced relative to each other con- 
nected by certain paths which are classified with respect 
to irreducible representations of the symmetry group of 
the hypercubic lattice. That approach is similar to the 
one adopted by Basak et al. [l3. 120L l2ll I22I] for baryons. 

In our paper we test derivative sources in a quenched 
excited meson spectroscopy calculation. In particular we 
study the pseudoscalar, scalar, vector and pseudovector 
channels. Depending on the channel we find considerable 
improvement of the signal for some of the excited and 
ground states. Preliminary results with our derivative 
sources were already reported in [23|, [24|, . 

II. SETTING OF THE CALCULATION 

A. Variational method 

The central idea of the variational method @, [1] is to 
use several different interpolators Oi,i — 1, . . . N with 
the quantum numbers of the desired state and to compute 
all cross correlators for interpolators projected to fixed 
spatial momentum (zero in this work), 

C(% - (Oi(t)O](0)>. (1) 

In Hilbert space these correlators have the decomposition 

C(t) tJ =J2(m\n){n\O]\0)e- tM ". (2) 

n 

Using the factorization of the amplitudes one can show Q 
that the eigenvalues X k (t) of the generalized eigenvalue 
problem 

C(t)v k = X k (t)C(t )v k , (3) 
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behave as 

A fe (t)oce- tM '=[l + 0(e- tAMfc )], (4) 
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where Mk is the mass of the fc-th state and AMk is the 
difference to neighboring states, fn Eq. the eigenvalue 
problem is normalized at a timeslice to <t. 

Equation (U) shows that each eigenvalue predomi- 
nantly decays with a single mass: The largest eigenvalue 
decays with the mass of the ground state, the second 
largest eigenvalue with the mass of the first excited state, 
and so on. Thus, the variational method disentangles 
the signals of the ground- and excited states. As a conse- 
quence simple, stable two-parameter fits become possible. 

At this point we remark, that the variational method 
also treats ghost contributions correctly, which in some 
channels show up in a quenched or partially quenched cal- 
culation at small quark masses (26| . It was shown in (2?| 
that in the variational approach the ghost contribution 
couples to an individual eigenvalue (up to the correction 
term) in the same way as a proper physical state. Thus, 
ghost contributions are disentangled from the physical 
states and need not be modeled in the further analysis of 
the exponential decay of the eigenvalues. 

To visualize the results and to determine possible fit 
ranges, we plot effective masses which are built from the 
ratios of eigenvalues 

•M t + 5)- ta (x£TIj)- <5) 

For those values of t where the exponential decay of the 
eigenvalue is governed by a single state, the effective 
masses form pronounced plateaus. 

It is an interesting observation that for the same values 
of i, where the effective mass plateaus form, also the 
corresponding eigenvectors are approximately constant 
as a function of t. An example of this behavior is given 
in Fig. O (discussed later), where we show the entries 
of the eigenvectors for the three largest eigenvalues as a 
function of t. This time-independence of the eigenvectors 
serves as a "fingerprint" for the physical states. To be 
more precise, the eigenvectors we use for such fingerprints 
are the eigenvectors of the regular eigenvalue problem 

C(t )-iC(t)C(t )-ivi = X k (t)vl , (6) 

which obviously has the same eigenvalues A&(t) as the 
problem ([3]), but gives rise to orthogonal eigenvectors v' k . 

In our analysis we only fit states which give rise to a 
plateau in both the effective mass and the corresponding 
eigenvector. Ideally, to should be chosen large. However, 
large to also tends to increase the statistical noise. We 
explored the dependence of the effective mass plateaus 
on the timeslice to and usually found the best results for 
to = 1. So unless noted otherwise, t Q will be fixed to 
to = 1 for all results quoted (our sources are located at 
t = 0; see next section). 

B. Smeared sources and sinks 

In order to optimize the overlap with the ground and 
first few excited states, one commonly uses quark smear- 



ing. We first construct extended sources by Jacobi smear- 
ing [2^,[2^| of point sources So located at timeslice t = 0: 

S { a ' a) (y,t) p , c = S(y,0)6(t,0)6 pa 6 ca , (7) 

N 

S (a,a) = K n H n S ( o a ' a) , (8) 
n=0 
3 

ff(x,y)=^(f/ i (x,Q)5(f+i,y) 

i=l 

+U i (x-i,0)^S(x-i,y)^ . (9) 

The smearing has two parameters k and N and leads to 
gauge covariant, approximately Gaussian shaped sources 
of different width. We use the same combinations of pa- 
rameters as in Q and we refer to our sources as "narrow" 
(S n ) and "wide" (S w ). 

Our derivative quark sources Wg i are constructed by 
applying a covariant derivative to the wide sources: 

Pi(x,y) = Ui(x, 0)5(x + i,y) - Ui(x - i, 0)'S(x — i,y) , 

(10) 

W di = PiS w . (11) 

These sources are then used in the construction of meson 
interpolators of definite quantum numbers. 



C. Meson interpolators 

Table U shows our interpolators for the different meson 
channels considered. In the first column the interpolators 
are numbered according to their structure rather than 
consecutively. Numbers 1-6 denote the Jacobi-smeared 
interpolators of while the interpolators 7-12 contain 
at least one derivative. Notice also that only the combi- 
nation of number and quantum numbers uniquely labels 
an interpolator. 

In some cases, an (anti-) symmetrization of the inter- 
polators is necessary to obtain the correct behavior un- 
der charge conjugation. Therefore, interpolators denoted 
as ugTd n / w in Table Q] should be read as ugTdn^ — 
u n / w Tdg i . We restrict ourselves to light, isovector (1 = 1) 
mesons with degenerate quark masses m u = rrid- 

All interpolators have been classified by their contin- 
uum quantum numbers I J PC , both for non- vanishing 
quark mass and in the chiral limit. As usual, P is the 
spatial parity, J is the total spin, and / the isospin. For 
a neutral qq system [46J the C-parity is related to the 
other quantum numbers in a standard way. To simplify 
the notation, we omit the spin and isospin projections in 
the notation. 

In addition to the quantum numbers listed in table 
U the lattice interpolators will also couple to continuum 
states with higher J due to the loss of rotational sym- 
metry 

(HE!. The lattice interpolators for PC mesons 
couple also to J > 4, but this does not influence our 
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FIG. 1: Diagonal entries of the correlation matrix as a function of t. The data are for bare quark mass am q = 0.04. The 
(lhs. plot) and 1 (rhs.) channels are shown. The numbers next to the correlators are according to Table H] 



conclusions since there are no observed resonances with 
J > 4 in the energy regime of interest. The lattice in- 
terpolators for \ PC mesons couple also to J > 3 and the 
issue is discussed in section UlI CI 

In the chiral limit, the different interpolators listed in 
Table Q] can be classified into representations of chiral 
SU(2) L x SU(2) R and U(l) A groups [13, EH, as well as 
with respect to their partial wave 2S+1 Lj decomposition 
[321 ]. Below we review these properties. 

We label with R the index of the chiral representa- 
tion; R = (0, 0), (1/2, l/2) a , (1/2, l/2) 6 , or (0, 1) © (1, 0). 
The chiral basis {R; IJ PC } is obviously consistent with 
Poincare invariance. Each of the interpolators in Table U 
has a fixed U(1)a transformation property. Namely, all 
those interpolators that belong to (0, 0) or (0, 1) © (1, 0) 
of SU(2)l x SU(2)r are scalars with respect to U(1)a, 
i.e., they transform into themselves under a U(1)a trans- 
formation. However, interpolators with opposite spa- 
tial parity and the same spin J and isospin / from 
the distinct (1/2, 1/2) Q and (1/2, l/2) representations of 
SU(2)l x SU(2)r transform into each other upon U 

The set of quantum numbers {R; IJ PC } uniquely fixes 
a partial wave content |/; 2S+1 Lj) of the quark-antiquark 
system in the center-of mass frame [32j . In particular, the 
different interpolators from Table U with quantum num- 
bers 1 ++ belonging to (1/2, 1/2)5 represent the |1; 3 Pq) 
partial wave in the qq system, irrespective of the num- 
ber of derivatives in the interpolator. Note however that 
there are some interpolators in the 1 ++ channel (num- 
ber seven and eight) which will not couple to the scalars 
in the chiral limit, since they belong to the (0, 1) © (1, 0) 
representation which requires J > 1. 

For the 1 h sector, there are two types of interpo- 
lators: interpolators (1-3, 9-10, 11) which transform as 
(l/2,l/2) n and represent the ll; 1 So) partial wave and 
time components of pseudovector interpolators which 
couple to pseudoscalars due to PCAC and which belong 



to the (0, 1) © (1,0) representation. 

In the 1 1 ++ sector, all interpolators transform as 
(0, 1) ffi (1, 0) and couple only to the |1; 3 Pi) partial wave. 

However, there are two kinds of interpolators with 
quantum numbers 1 1 . They are the fixed and or- 
thogonal superpositions of two different partial waves: 

|(0,1) + (1,0); 11") = y||l;3 (Sl ) + y||l;3 i j l)i 
|(l/2,l/% 11") = yi|l;^ l) _y||l ; 3 Cl) . 

The interpolators 1-3 and 11 from Table [J belong to 
the |(0, 1) + (1, 0); 1 1 ) representation, while all others 
transform as |(l/2, l/2) 6 ; 1 1 — ). 

D. Technicalities 

For our analysis we used 99 uncorrelated quenched 
gauge configurations generated with the Liischer-Weisz 
gauge action [33,|34j]. We work on a 16 3 x 32 lattice with 
a = 0.148 fm determined [35| from the Sommer parame- 
ter (using ro = 0.5 fm). For comparison we also use old 
data from a 20 3 x 32 lattice at a = 0.119 fm, where, how- 
ever, only the Jacobi-smeared sources without additional 
derivatives are available. The boundary conditions for 
the gauge fields are periodic in all four directions. The 
quark propagators were computed from the Chirally Im- 
proved Dirac operator [36l \37\ with periodic boundary 
conditions in space and antiperiodic boundary conditions 
in the time coordinate. We study several quark mass 
parameters in the range am q = 0.02 . . . 0.2. We fold in- 
dividual entries of the correlation matrix resulting from 
propagation in positive and negative time direction ac- 
cording to their symmetry, which reduces the statistical 
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TABLE I: List of our meson interpolators. The numbers in the first column together with the quantum numbers I J given 
in the third column label the interpolators uniquely. The fourth and fifth column specify the chiral representation and the 
coupling in the chiral limit. The lattice interpolators may also couple to states with higher angular momentum J [T3 . [Isj ]. 
The subscripts n and w refer to narrow and wide smearing of u and d quarks. The subscript di denotes derivative smearing 
in i-direction. Where it appears, the index i is summed over the spatial directions 1,2,3. The time direction is 4 and the 
corresponding Dirac matrix is 74. For the vector and pseudovector channels the index k (not summed !) can have values 
fc = l,2,3. 
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errors and improves the quality of the data significantly. 
Unless noted otherwise, the errors we quote are statistical 
errors determined with the Jackknife method. 

Where possible we also indicate the systematical un- 
certainties by a shaded band (Figs.[3l[4]etc.). The upper 
and lower limits of this band are obtained by repeating 
the fits of the eigenvalues using different fit ranges and 
varying the interpolators used in the correlation matrix. 
Although this is certainly only a rough estimate of the 
systematic uncertainty in the mass determination, we re- 
fer to these error estimates as "systematic errors". We 
stress, however, that these errors do not include the error 
introduced by the quenched approximation. 

For plots with the pion mass squared on the horizontal 
axis, we use a specific combination of gaussian interpo- 
lators for extracting the ground state pion mass. The 
corresponding statistical error gives rise to the horizon- 
tal error bars in some of our plots which are, however, 
smaller than the symbols used. 

When fitting eigenvalues obtained with the variational 
method we experimented with both correlated and un- 
corrected (two parameter) exponential fits. For the cor- 
related fits we used a jackknife estimate of the correlation 
matrix which was not always stable with our ensemble of 
configurations. For the generation of the plots we there- 
fore resorted to simple uncorrelated fits throughout, us- 
ing only the diagonal elements of the covariance matrix. 
For the cases, where correlated fits are stable and a direct 
comparison with uncorrelated fits is possible we find that 
the latter give larger statistical errors [38| . This explains 
the rather small x 2 /d.o.f. we find. Thus the uncorre- 
lated errors we give in the plots and in the appendix are 
probably overestimated. 



E. A first look at the derivative sources 

Figure[T]shows the diagonal elements of the correlation 
matrix for some of the channels considered. Compared 
to the interpolators from reference Q (without derivative 
sources), the interpolators with derivative sources show 
stronger contributions from excited states, i.e., they have 
a steeper slope for small Euclidean time t. Nevertheless, 
for all interpolators the ground state in the respective 
channel dominates the behavior at large time separation, 
i.e., all correlators in Fig. Q] display the same slope at 
sufficiently large t. 

We remark that the propagators in the ++ channel 
at the lightest quark masses show a deviation from that 
pattern for small t due to ghosts. A more detailed dis- 
cussion of the scalar correlators will be given in Section 

mm 

As suggested by Fig. [TJ the pion ground state can 
be fitted from a single diagonal correlator at all quark 
masses. Stable plateaus are obtained in the time interval 
t = 6 ... 15 where cosh- fits can be performed for the indi- 
vidual correlators. For the lowest quark mass am q = 0.02 
the statistical error usually is of the order of 1 — 2% of the 
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FIG. 2: First three eigenvalues from the generalized eigen- 
problem ([3} for the channel as a function of t. 

fitted value with the derivative sources leading to a some- 
what larger error. An exception is interpolator 12 where 
the statistical error is about 3.5% of the fitted value. At 
larger masses the error is substantially smaller. The fit 
values for the ground state pion mass from the different 
individual correlators, shown in the lhs. plot of Fig. [TJ 
agree within two sigma. 

For the 1 _ ~ channel similar observations hold. 



F. Contributions from backward propagation 

In Fig. [5] we show the three largest eigenvalues obtained 
from the generalized problem {3| for the pion case. It 
is obvious that on the logarithmic scale used in this plot 
the three eigenvalues give rise to three essentially straight 
lines for sufficiently small t. The different slopes corre- 
spond to the masses of the ground-, first- and second 
excited states. 

For mesons, forward and backward propagation be- 
haves in the same way. Each interpolator coupling to 
a particular state at early times will also couple to the 
same, but backward running, state at later time. Higher 
channels in the generalized eigenvalue problem may in 
this way turn into lighter state signals. This interesting 
observation (see also [3!|), particular for the generalized 
eigenvalue problem, can be made for the second eigen- 
value in Fig. [5J At t ~ 9 (when using to = 1) the data 
points of the second eigenvalue change their behavior and 
start to increase again. Beyond t = 9 the data points 
form a straight line with positive slope. This upward 
pointing straight line turns around again at t ~ 13 and 
from there on decays with a slope corresponding to the 
ground state mass. Also the slope of the upwards point- 
ing piece between t = 9 and i = 13 (which is then con- 
tinued by an upward pointing piece of the largest eigen- 
value), has the slope of the ground state mass. 

This avoided level crossing scenario has an important 
consequence: The generalized eigenvalue problem disen- 
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FIG. 3: First excited state of the pseudoscalars for two dif- 
ferent sets of interpolators. The error bars are statistical only 
and the shaded region indicates the additional systematic er- 
ror as discussed. The filled symbol corresponds to the exper- 
imentally measured 7r(1300). 
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FIG. 4: Second excited state for the pseudoscalars. The filled 
symbol indicates the 7r(1800). 



tangles the forward propagating ground and excited state 
masses only up to the first crossing with the backward 
running lightest propagator, which for our example hap- 
pens at t ~ 9. Beyond that value also the second eigen- 
value, which for small t is dominated by the first excita- 
tion, couples to the lighter ground state (running "back- 
wards" ) and no longer provides information on the exci- 
tations. In particular in the h channel, where a back- 
ward running light pion crosses with the second eigen- 
value already at a small i, this effect limits the analysis 
of excited states [43] . 

The comments of this subsection are particularly im- 
portant for light pseudoscalars. For the pions it is the 
backward running contributions that limit the fit range 
with the generalized eigenproblem, leading to errors com- 
parable with the errors obtained from fitting single cor- 
relators with the correct functional form. 



III. RESULTS FOR INDIVIDUAL MESON 
CHANNELS 

A. The CT+ channel 

Figure [3] shows the results for the first excited state 
of the pion. We display the results for two different sets 
of interpolators. Circles are used for the combination of 
gaussian interpolators 1, 4, 5, 6, while the squares corre- 
spond to the combination 1, 4, 6, 9, 12 of gaussian and 
derivative interpolators. While combinations of gaus- 
sian and derivative interpolators allow for fits at slightly 
lighter quark masses, the interpolators with derivative 
sources couple weaker to the ground state. The system- 
atical uncertainty (shaded region in the figure) from the 
choice of interpolators is consistent with statistical effects 
as reasonable combinations of four or more interpolators 
lie within one sigma of our final fit result. While there is 
no significant improvement, the new results with a larger 
basis nicely confirm the existence of the measured state. 

With a combination of gaussian and derivative inter- 
polators, it is also possible to obtain fits for a second 
excited state which could not be observed before. This 
state is displayed in Fig. 2J In the chiral limit, this state 
can most likely be identified with the tt(1800). Fits with 
various different combinations of interpolators lead to the 
same results which all show stable eigenvector entries. 

It is instructive to look at the components of the eigen- 
vectors for all three states observed in the pseudoscalar 
channel. Figure [5] shows such a plot. While the deriva- 
tive interpolators 9 and 10 do not contribute significantly 
to the ground and first excited states, they are most im- 
portant for obtaining the newly observed second excited 
state. This behavior is qualitatively the same for all pos- 
sible combinations of interpolators where the second ex- 
cited state could be seen. 

While the reduction in the statistical error for the first 
excited state can merely be attributed to an enlargement 
of the basis, the second excited state is only observed 
when including derivative interpolators. We would like to 
stress that a correlation matrix of similar size consisting 
solely of non-derivative operators does not enable us to 
see this excitation. 



B. The ++ channel 

In the 1 ++ channel contributions from ghosts [26|, 
arising from the 77/77; contribution to 
the isovector-scalar correlators, are expected and must 
be identified for a clean interpretation of the data. These 
unphysical contributions due to quenching have a nega- 
tive spectral weight and dominate the correlators at small 
quark masses, leading to correlators which become nega- 
tive at intermediate time separations. Figure [6] shows the 
diagonal correlators for the scalar channel at our smallest 
quark mass am q — 0.02. While some of them (1-3) dis- 
play a very prominent ghost contribution, others show 
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FIG. 5: Eigenvector components from the standard eigenvalue problem |6| as a function of t. The eigenvectors correspond 
to the pion ground state (lhs. plot), the first and second excited states (center and rhs. plots) and are obtained from a 4 x 4 
matrix of the interpolators 1, 6, 9, 10, at quark mass am, = 0.12. 



a much smaller contribution. Correlators 7 and 8 feel 
no effects from ghosts at all with our limited statistics, 
which may be related to their different chiral structure 
(see Table QJ . We note that these two interpolators do 
not seem to couple in the dynamical case, as our prelim- 
inary dynamical results do not show a signal for these 
correlators. 

While the variational approach enables us to disentan- 
gle most of these ghost contributions [271] . at low quark 



masses the quality of the observed plateaus quickly de- 
teriorates for the gaussian interpolators leading to large 
error bars for the ao ground state. 

The additional interpolators with derivative sources 
enlarge the correlation matrix, which is vital in the pres- 
ence of ghost contributions. Furthermore, it is these in- 
terpolators which couple only weakly to the ghosts. Us- 
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FIG. 7: Ground state mass of the isovector, scalar (ao). The 
error bars are statistical only and the shaded region indicates 
the additional systematic error described in the text. The 
black circle shows the (experimentally) measured ao(1450) 
and the triangle indicates the ao(980). 



ing the variational method we are able to disentangle the 
leading ghost contribution and in some cases even a sub- 
leading ghost contribution, corresponding to a rj' n state 
with relative momentum. 

Figure [7] shows the results for the largest eigenvalue 
of the variational analysis. The plot demonstrates that 
derivative sources enable us to perform fits at smaller 
quark masses with reduced statistical errors in the inter- 
mediate and heavy quark mass region. 

The shaded region in Fig. [7] indicates our estimate 
of the additional systematic errors due to the choice of 
fit ranges and interpolators considered. As it remains 
unclear which systematics cause the dependence on the 
choice of interpolators in this channel, we refrain from an 
extrapolation to the physical mass region. Interpolators 
with strong ghost contributions, however, tend to lead to 
higher mass values suggesting the ground state to be in 
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the region of the ao(1450), while a fit with those inter- 
polators containing no visible ghost contribution leads to 
values in the lower parts of the shaded region. 

For this channel large quenching effects might influence 
the result. Thus it is not clear whether the data should 
extrapolate to the ao(980) or the ao(1450). More light 
will be shed on this channel only when dynamical data 
are available where no quenched ghosts are present [44j . 
If the systematic deviations we observe in the choice of 
interpolators are due to ghosts, one should expect the 
results to become more consistent with dynamical data. 

A review of issues faced in the scalar channel and a 
detailed discussion of the possible nature of the isovector 
scalar groundstate can be found in the recent review by 
McNeile 



C. The 1 channel 

As mentioned before, the ground state for the vector 
meson, the p(770), can also be fit from single correla- 
tors. However the results improve quite drastically when 
a matrix of interpolators is used. 

Figure [8] shows the 1 meson ground state and il- 
lustrates the good quality of the data in the quenched 
approximation, where no decay is possible. The results 
from different interpolators and fit ranges agree within 
error bars. 

The interpretation for the first and second excitation 
in the 1 channel is less clear. From experiment we 
know of multiple excitations with J = 1 below 2 GeV, 
the most established being the p(1450) and the p(1700). 
In addition, due to loss of continuous Lorentz symmetry, 
some of the lattice interpolators we chose may in prin- 
ciple couple to continuum states with higher J [Til [l8| , 
and there is at least one such excitation known in the 
vector channel, the p3(1690). Excluding this possibility 
would be a difficult task which might be overcome by 
taking a look at different irreducible representations of 
the hyper-cubic group where degeneracy of states in dif- 
ferent representations can be used to identify the correct 
J, as has been demonstrated for baryons fig . I20I l2l| . 

As can be seen in Figs. [9] and [TOl the values obtained 
from combinations of gaussian interpolators agree qual- 
itatively with the values obtained from the larger basis, 
while the larger basis leads to overall smaller error bars 
and somewhat more stable plateaus. Nevertheless, the 
two excitations are too close together and would both be 
consistent with the o(1700). This problem has already 
been encountered in [9j, where both the 16 3 x 32 lattice 
with spacing a = 0.148 fm, and a finer 20 3 x 32 lattice 
with a = 0.119 fm have been used. The data from the fine 
lattice lead to two distinct excitations compatible with an 
interpretation as the physical p(1450) and p(1700). To 
demonstrate this, we also included an alternative fit of 
the old data from the fine lattice in Figs. [91 and [TOl 

Inspecting the eigenvectors of the states, we can iden- 
tify them by their operator content and come to the con- 
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FIG. 8: Ground state for the p(770) meson. 
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FIG. 9: First excited state for the p meson, compared to the 
first experimental meson resonance with J = 1, the p(1450). 
The diamonds represent an alternative fit of data from [j| for 
a finer lattice of the same volume (20 3 x 32, a = 0.119 fm). 
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FIG. 10: Second excitation in the vector channel compared 
to the p(1700). Again the diamonds indicate data from the 
finer lattice. 
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elusion that they are indeed the same on both lattices. 
Moreover, as both lattices have the same physical vol- 
ume, this leads us to an interpretation of the difference 
as a discretization effect. Such an explanation seems rea- 
sonable as the states are rather close to each other. We 
therefore are confident that the method works and that 
the same calculation on finer lattices would lead to results 
in better agreement with experiment. 



D. The 1 H 



channel 



For the 1 ++ channel, there is only one interpolator con- 
taining derivative quark sources. Figure [TT] shows data 
obtained from different combinations of gaussian interpo- 
lators and the one containing derivatives. An indication 
of error bands, as shown for the other channels, has been 
omitted here, since the two combinations plotted already 
show the extremes. 

Figure [TT] demonstrates a clear improvement in the de- 
scription of the ground state using the interpolator with 
derivative quark sources. While the results from the 
gaussian and the full sets agree qualitatively, the statisti- 
cal errors towards smaller quark masses are significantly 
reduced. The reason for this are longer, more stable ef- 
fective mass plateaus allowing for larger fit ranges. At 
larger quark masses there is a slight deviation of the or- 
der of two sigma. 

The excited state previously observed stays the same if 
one includes the new interpolator in the analysis. Look- 
ing at the components of the modified eigenvalue problem 
we see that this interpolator contributes only weakly to 
this excited state. 
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IV. SUMMARY 

In this article we have explored the impact of deriva- 
tive sources in light meson spectroscopy. The sources are 
obtained by applying a covariant derivative on a Jacobi- 
smeared quark source. Interpolators based on derivative 
sources were tested in a quenched excited meson spec- 
troscopy calculation based on the variational method. 

We find that both, ground- and excited state signals 
may be improved, depending on the channel. For the 
h channel we find that the quality of the first excited 
state improves and fits become possible for smaller quark 
masses. In addition a second excited state can be iden- 
tified when adding interpolators with derivative sources 
in the correlation matrix. 

The ++ channel is dominated by the presence of 
ghosts at small quark masses which we essentially disen- 
tangle using the variational method. The additional in- 
terpolators with derivative sources are helpful since they 
enlarge the correlation matrix and some of them couple 
only very weakly to the ghost states. They allow one 
to fit scalar masses at lower quark masses, although the 
results still depend significantly on the choice of inter- 
polators used in the correlation matrix. This indicates 
that results for this channel from a quenched calculation 
should be interpreted only with the necessary caution. 

For the 1 mesons we demonstrate that the inclu- 
sion of derivative sources leads to results with a smaller 
statistical error which agree excellently with the results 
published before. 

For the 1 ++ channel we show that the derivative 
sources drastically improve the signal for the ground 
state. At the same time a matrix of interpolators en- 
ables us to identify a second excited state. 
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interpolators 


m q 


mass [GeV] 


error 


fitrange 


X 2 /d.o.f. 


1,4,5,6 


0.20 


1.0535 
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4-10 


0.09 


0.16 
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0.0025 


4-10 


0.11 


0.12 


0.8098 
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4-10 


0.13 


0.10 


0.7391 
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0.0029 
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0.0029 
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0.04 
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0.03 
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4-10 


0.12 


0.02 


0.3420 
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4-10 


0.13 


1,6,9,10 


0.20 


1.0536 
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3-9 
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0.9380 
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0.12 
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0.06 
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0.05 


0.5274 
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0.03 
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1,4,6,9,12 
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0.0024 
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0.12 


0.8096 
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3-9 


0.17 


0.10 
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0.0028 
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interpolators 


m q 


mass [GeV] 


error 


fitrange 


X 2 /d.o.f. 


1,4,5,6 


0.20 


1.974 


0.062 


3-7 


0.10 


0.16 


1.891 


0.069 


3-7 


0.09 


0.12 


1.813 


0.079 


3-7 


0.04 


0.10 


1.778 


0.086 


3-7 


0.02 


0.08 
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0.088 


3-6 


0.01 


0.06 


1.698 


0.100 


3-6 


0.02 


1,6,9,10 


0.20 


1.999 


0.047 


3-7 


0.02 


0.16 


1.917 
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3-7 


0.01 


0.12 


1.834 


0.062 


3-6 


0.07 


0.10 
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0.14 


0.08 
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0.21 


0.06 
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0.05 
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0.05 


1,4,6,9,12 


0.20 


1.972 


0.064 


3-7 


0.11 


0.16 


1.888 


0.069 


3-7 


0.10 


0.12 


1.805 


0.078 


3-7 


0.03 


0.10 


1.766 


0.084 


3-7 


0.05 


0.08 


1.730 


0.088 


3-6 


0.06 


0.06 
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0.104 


3-6 


0.03 


0.05 


1.646 


0.117 


3-6 


0.07 


0.04 
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3-6 


0.04 


TABLE III: First excited state in the + channel. 
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m q 


mass [GeV] 


error 


fitrange 


X 2 /d.o.f. 


1,6,9,10 


0.20 
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0.01 


0.16 
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0.02 
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0.10 
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interpolators 


m q 


mass ]GeVJ 
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fitrange 
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0.05 
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11 
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0.08 


1.554 
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TABLE IX: Ground state in the 1 ++ channel. 
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0.10 


2.241 
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